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Abstract 

A previous stability condition (see Throumoulopoulos and Tasso, 
Physics of Plasmas 14, 122104 (2007)) for incompressible plasmas with 
field aligned flows is extended to gravitating plasmas including self- 
gravitation. It turns out that the stability condition is affected by 
gravitation through the equilibrium values only. 
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The system under consideration is a constant density plasma with an 
incompressible flow aligned with the magnetic field, and maintained in a 
magnetohydrodynamic equilibrium by the Lorenz force and external and self- 
gravitation. It is governed by the following equations 

V = AB (1) 

with A an arbitrary function, 

V ■ V = 0, (2) 

V ■ Vv = J X B - VP + V0, (3) 

where v is the velocity, B is the magnetic field and J = V x B is the current 
density, (f) is the gravitational potential. The equations are written in con- 
venient units and the constant density is set to 1. (In fact, as will be proved 
below it suffices to assume constant density only at equilibrium.) In addi- 
tion we assume the existence of well defined equilibrium magnetic surfaces 
labelled by a smooth function %p. These equations differ from the correspond- 
ing equations of Ref. [Ij by the inclusion of the gravitational potential. 

Equations (1) and (2) lead to A = A('?/'), and using the identity v ■ Vv = 
- V X V X V, (l)-(3) lead to 

(1 - A^) J X B = V(P + ^ - 0) - P^Vy , (4) 

where B and v are the magnetic field and velocity field moduli respectively. 
Taking the scalar product of Eq.(4) with B gives 

P + ^-(^ = PsM^ (5) 

where Ps^{ip) is the arbitrary function in the presence of gravity. 
Consequently, Eq.(4) can be put in the form 

(1 - A^) J X B = P;V^ - (A2)'^VV^, (6) 

with 

JxB = (?^(^,p2)vv., (7) 



where 
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The only difference with Ref. pQ is the presence of the gravitational po- 
tential in (4)- (8). This potential obeys Poisson's equation 



A0 = -1, (9) 

whose solution is the superposition of the internal potential 0j given by 
(^i = Jy (x-x')l '^^'^^' where V is the volume of the plasma considered, and 
an external potential 0e due to fixed masses outside the plasma. 

Let us now come to the stability of our flowing and gravitating equilibria. 
The general problem without gravity has been already formulated in Ref.[2] 
in terms of the Lagrangian displacement ^. Its application to equilibria with 
field aligned incompressible flows and homogeneous density was investigated 
in Ref . [Ij and Ref . [3J . The correct derivation of the final condition was done 
in Ref.jlJ. This condition states that a general steady state of a plasma of 
uniform density and incompressible flow parallel to B is stable to small three- 
dimensional perturbations if the flow is sub-Alfvenic (A^ < 1) and A > 0, 
where 

A = -/[(l-A2)(Jx V^)-(B-V)V^+^^|V^r(V^-V^+^?|VV^r)]. (10) 

The quantity g is like the of Eq.(8) but the Pg in it does not contain the 
gravitational potential 0. Now, A > remains valid if g in (10) is replaced 
by g^i^. The reason is that the perturbed gravity in the linearized dynamical 
equations vanishes for the V ■ ^ = and homogeneous density at equilibrium. 
To prove this statement we employ the relations[2] 

Pi = -V ■ (p^) , (11) 

vi = § + (v-V^-^-Vv), (12) 

where the subscript 1 denotes first order quantities. Acting the divergence 
operator on f|T2|) and using the identity V x (^ x v) = v ■ — ^ ■ Vv + 
^ (V ■ v) - V (V ■ ^) with V ■ V = yields 

V-vi = V-^ + vV(V-^)^^(V-0. (13) 

Therefore, if V ■ ^ = it follows V ■ vi = and, inversely, V ■ vi = implies 
V ■ ^ = for t > to if this relation is satisfied initially at t = to- Then, ( II ip 
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in conjunction with the hnearized continuity equation, dpi/dt + v • Vpi = 0, 
lead to pi = 0. The fact that for uniform equihbrium density, the density 
remains uniform in the perturbed state because of incompressibihty was not 
noticed in Refs.jl] and p!]. 

We conclude by stating that the above mentioned condition is sufficient 
for stability if we replace g in (10) by [Eq. (8)], which means that this 
stability condition is affected by gravity through the equilibrium values only. 
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